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ABSTRACT 

The stability properties of a low density ultra relativistic pair beam produced in the intergalactic 
medium by multi-TeV gamma-ray photons from blazars are analyzed. The problem is relevant for 
probes of magnetic field in cosmic voids through gamma-ray observations. In addition, dissipation of 
such beams could affect considerably the thermal history of the intergalactic medium and structure 
formation. We use a Monte Carlo method to quantify the shower properties, in particular the bulk 
Lorentz factor and the angular spread of the beam, as a function of distance from the blazar. We 
find that the fastest growing modes, like any perturbation mode with even a very modest compo- 
nent perpendicular to the beam direction, requires a kinetic treatment. Combined with the effect of 
non-linear Landau damping, which suppresses the growth of plasma oscillations, the beam relaxation 
timescale is found typically longer than a Hubble time. Finally, density inhomogeneities associated 
with cosmic structure induce severe loss of resonance between the beam particles and plasma oscilla- 
tions, strongly inhibiting their growth. We conclude that relativistic pair beams produced by blazars 
in the intergalactic medium are stable on timescales much longer compared to the electromagnetic 
cascade's. There appears to be no effect of pair-beams on the intergalactic medium. 
Subject headings: gamma rays: general - instabilities - intergalactic medium - plasmas - radiation 
mechanisms: non-thermal - relativistic processes 



t— 1 
00 

o 



X 



1. INTRODUCTION 

Streaming relativistic particles are common in tenuous 
astrophysical plasma and their propagation and stability 
properties a recurrent theme. Examples inc lude transien t 
hard X-ray e mission in the solar corona (lBen7|[l99l . 
quasars' jets (jLesch fe Schlickeiser 1 119871). cosmic-rays 
strea ming out of star forming galaxies (jMiniati fc Bell I 
Hoiil) and cosmic-ray trans port in the intracluster 
medium (|Ensslin et al. I [201 If) . Propagation and stabil- 
ity properties, related in particular to the exitation of 
plasma waves, is subject of attentive investigation as they 
can play a crucial role in the interpretation of observa- 
tional data. 

Ultra-relativistic beams of e ± pairs are also generated 
in the intergalactic medium (IGM) by very high energy 
gamma-rays from distant blazars, by way of photon- 
photon inte ractions with the extragalactic background 
light (EBL, ICould fc Schreder lfT967t ISchlickeiser elaTI 
20121). While blazars' spectra, and in particular their 
multi-TeV cut-off features, have been studied in d etail 
to constraint the EBL (e.g. lAharonian et al. |[2006l ). re- 
cently multi-GeV and TeV blazars observations have also 
been used to c onstrain magnetic field in cosmic voids for 
the first time (|Neronoy fc Vovk |[2010t iTavecchio et al. I 
120101 iDermer et al. 1 120111 : iTavlor et al. II2011D . In fact, 
for flat enough blazar's spectra, the electromagnetic cas- 
cade should produce an observable spectral bump at 
multi-GeV energies. The absence of such a feature in 
a number of observed blazars is ascribed to the presence 
of a sufficiently strong magnetic field, B v > 10 _18 G, to 
deflect the pairs in less then an inverse Compton length, 



lie - Mpc (E+ ZTeVy Ul + z)~ 4 , where z is the cosmo - 
logical redshift ( Pla ga|[T995t IN eronov fc Semikoz ll2009h . 

However, in principle the pair-beam is subject to 
various instabilities, in particular microscopic plasma 
instab ilities of the two-strea m family. On this ac- 
count, iBroderick et al. I (|2012D conclude that transverse 
modes of the two-stream instability act on much shorter 
timescales than inverse Compton scattering, effectively 
inhibiting the cascade and invalidating the above mag- 
netic field measurements. In addition, as a result of 
the beam's relaxation, substantial amount of energy 
would be deposited into the IGM, with dramatic con- 
sequences for its ther mal history (jChang et al. 1 120121 : 
iPfrommer et al. 112012ft . 

In this paper, we reanalyze the stability of blazars in- 
duced ultra-relativistic pair beams. In particular, we 
use a Monte Carlo model of the electromagnetic shower 
to quantify the beam properties at various distances 
from the blazar, and analyze the stability of the pro- 
duced beam following the w ork of Breizman, Rytov and 
collaborators (r eviewed in iBreizman fc Rvutov I 119741 : 
lBreizman~lll990ft . We pay particular attention to the role 
of small, but finite, angular spread of the beam. We find 
that even for very modest perpendicular components of 
the wave-vector, the analysis of the instability requires 
a kinetic treatment. In addition, when non-linear Lan- 
dau damping is taken into account, the growth rate of 
plasma oscillations is suppressed, so that the beam re- 
laxation timescale becomes longer than a Hubble time. 
Finally, the resonance condition cannot be maintained in 
the presence of density inhomogeneities associated to cos- 
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mological structure formation, which also act to dramat- 
ically suppress the beam instability. Non-linear Landau 
damping was not explicitly included in Broderick's calcu- 
lations and density inhomogeneities not considered, but 
they both appear to considerably stabilize the beam dur- 
ing a Hubble time, thereby ruling out effects of blazars' 
beam on the thermal history of the IGM and structure 
formation. 

The rest of this paper is organized as follows. Sec. [2] 
summarizes the physical properties of pair beams pro- 
duced by blazars and present the results of the Monte 
Carlo model. The two-stream instability in both the re- 
active and kinetic regimes is discussed in Sec. [3J Nonlin- 
ear effects are discussed in Sec. [4] where the timescales 
for the relaxation of the pair beam are also derived. Fi- 
nally, Sec. [5] briefly summarizes the results. 

2. PAIR BEAMS IN VOIDS 

The unstable behavior of an ultra-relativistic pair 
beam depends, amongst others, on its density contrast 
to the IGM, Lorentz factor, angular and velocity spread. 
These characteristic quantities derive from the energy 
and number density of the pair producing photons, i.e. 
the blazar's spectral flux, F 7 , and the EBL model 
The proc ess has been studied extensively in the litera- 



ture ( e.g. iGould &: Schreder I 19671 : iBonometto fc Reesl 
119711 : Schlickc iser et al. 112012 ) and in the following we 



briefly summarize its qualitative features before we de- 
scribe the results of our Monte Carlo model of a blazar 
induced cascade. 

2.1. Basic Qualitative Features 

Pairs arc most efficiently created just above the energy 
threshold for production, i.e. where 



s = £ 7 £ E bl(1 - cos<p)/2m 2 e c 4 > 1, 



(1) 



with <j) the angle between the interacting photons, E 1 and 
-Bebl the energy of the incident and target EBL photon, 
respectively, and the relativistic invariant, s, the center 
of mass energy square in units m 2 c A . The mean free 
path fo r the process depends on the details of the EBL 
model (jKneiske et al. 1 12004 iFranceschini et al. 1 120081) 
but is approximately 
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(l + z)- c Gpc, 



(2) 



with C = 4.5 for z < 1, C = other - 

wise (N cronov fc Semikoz 1 12009: Brodcric k et al. If2012h . 
The particle number density of the beam, rib, is set by 
the balance of pair production rate, 2F 1 jl 111 evaluated 
close to production threshold, and energy loss rate. If 
inverse Compton losses dominate then, at a distance D 
from the blazar such that F 1 = L^/AttD 2 , 
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where, i? 7 L 7 , is an estimate of the blazar's equivalent 
isotrop ic gamma-ray luminosi ty for a source at distance 
D (see. lBroderick et al. 1120121) . Each pair particle carries 
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Fig. 1. — Top: Peak (dash) and median (solid) pair beam energy 
as a function of distance from the blazar. Vertical bars correspond 
to 68% percentile of the energy spread about the median. A power 



law, dn-y/dEj oc , in the range 10 3 < E-y/r, 



< 10 8 , 



sumed for t he blazar's spectru m, together with the nominal EBL 
model from Aharonian (2001). Bottom: Angular distribution of 
beam pairs generated at 55 Mpc (solid), 350 Mpc (dash) and 1 Gpc 
(dot-dash) from the blazar. The colored curves indicate the con- 
tribution to the angular distribution for 350 Mpc from pairs in the 
energy range: 10 4 -10 5 (red), 10 5 -10 6 (yellow), 10 6 -10 7 (green), 
10 7 — 10 8 (blue), respectively. 

about half the energy of the incident gamma-ray, so the 
beam Lorentz factor is 



r, 



E-, 



beam 



2m P c 2 
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(4) 



Another important characteristic quantity of the beam is 
its angular spread, A#fc eam , determined by the distribu- 
tion of angles 9 between of pair produced particles with 
respect to the parent photon direction. This can be found 
to be related to Tb eam and the relativistic invariant as 
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2.2. Monte Carlo Model 



(5) 



The characteristic quantities of a pair beam gener- 
ated by gamma-rays from a blazar were computed us- 
ing a Monte Carlo model of the electromag netic cascade. 
The M onte Carlo code is fully described in lElviv et al. I 
(2009). The calculation assumes the nominal model of 
EBL in the range from 0.1 to 1000 fim (jAharonian 1120011) 
and a typical gamma-ray power law photon distribution, 
dn 7 /dE 7 oc E~ 2 , in the range 10 3 < E 1 /m e c 2 < 10 s , for 
the blazar's spectrum. 

Results for the pair beam generated according to the 
Monte Carlo model of the electromagnetic cascade are 
illustrate in Fig.Q] The top panel shows the peak (dash) 
and bulk (solid) energy of the pair beam, in units of 
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m e c 2 , as a function of distance from the blazar. The 
typical Lorentz factor of the pairs is a monotonically de- 
creasing function of distance, in the range r=(5-l)xl0 5 
up to a Gpc away from the blazar. This shows that 
pair production typically peaks at near infrared-optical 
{Eebl — 0.1-1 eV) EBL target photons interacting with 
gamma-rays with energy E 1 ~ (O.l-l)TeV. The verti- 
cal bars indicate the energy range encompassing 68% of 
the particles. In fact, there is considerable energy spread 
about the mean value, which decreases towards larger 
distances as the energy range of gamma-ray photons in- 
teracting with the EBL is also reduced. However, the 
beam particles remain always ultra-relativistic, a detail 
relevant in the analysis below. 

The bottom panel shows the angular distribution of 
beam pairs at distances of 55 Mpc (solid), 350 Mpc (black 
dash) and 1 Gpc (dot-dash) from the blazar. The angu- 
lar distributions dn/dO were simula t ed us ing analytical 
expressions from ISchlickcisc r et al. I (|2012t ) for angles in 
double photon collision. The beam angular spread is of 
order A8 ~ 10~ 5 , consistent with Eq. ([5]) and the value 
of r estimated above. The colored, dash curves show the 
angular distribution of pairs at 350 Mpc (dash) in four 
different energy bins, indicating that the beam angular 
spread is primarily determined by low energy pairs with 
T ~ 10 4 -10 6 (red and yellow). 

2.3. IGM in Voids 

Free electrons in cosmic voids are characterized by a 
number density n v ~ 2 x 10~ 7 (l + <5)(l + z) 3 cm -3 , where 
the typical over-den sity is S < —0.1, and a temperature 
T v ~ a few xlO 3 K (IPeebleslfl993l) . 

3. BEAM INSTABILITY: REACTIVE VS KINETIC 

The blazar induced pair beam is subject to microscopic 
instabilities, in particular two-stream like instabilities, 
of both electrostatic and electromagnetic nature. The 
beam is neutrally charged, so a return current is not 
induced. In the absence of magnetic field, luh <C oj p , 
with ojh the cyclotron frequency, u> p = (47rn v e 2 /m e ) 1 / 2 
the plasma frequency of the IGM in voids and e the 
positron's charge, the instability is predominantly asso- 
ciated to Cherenkov emission of Langmuir waves, which 
operates under the following resonant condition 



0. 



(6) 



where k is the wave- vector of the perturbation mode and 
v the beam particles velocity. The pair particles con- 
tribute equally to the dielectric function, as they have the 
same mass, number density, velocity distribution, and 
plasma frequency, uj p ^ = (47rn( ) e 2 /m e ) 1 / 2 . After sep- 
arating the contributions from the background plasma 
and the beam particles, the dispersion relation for Lang- 
muir waves_j_valid in the relativistic case, can be written 
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dp = 0. 



(7) 



where, /(p), is the distribution function of the beam par- 
ticles. There are two important regimes that character- 
ize the unstable behavior of the beam, namely reactive 
and kinetic. In the reactive case, the beam's velocity 
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Fig. 2. — Normalized growth rate, 74, /ttu p (nj, /n v ) , from Eq. 
in the plane fcyc/ojp — k±c/uj p . Bright is positive, dark is negative 
and the uniformly colored region is zero, as it lies outside the res- 
onant region. 

spread, Av, is negligible so all particles can participate 
to the unstable behavior and the growth rate of the in- 
stability is therefore fastest. In the kinetic regime, on 
the other hand, the velocity spread is considerable and 
only the resonant particles contribute to the growth of 
Langmuir waves, so the growth rate is slower than in the 
reacti ve case. Formally, the reacti ve regime is applicable 
when (|Breizman fc Rvutov II197IT ) 



|k • Av| < 7,. 



(8) 




where j r is the reactive growth rate. In this case, the 
integral in Eq. ([7]) can be solve in a simplified way, which 
involves neglect of the velocity spread around the mean 
value. This leads to the estimate of the reactive growth 
rate which, maximized al ong the longitudinal c omponent 
of the wave- vector reads (Fainbcrg et al.lll970T) 



(9) 



with fcii = uip/v, and ku, k± the components of the 
wave-vector parallel and perpendicular to the beam di- 
rection, respectively. It is well known that, for an ultra- 
relativistic beam (r 3> 1), the fastest growing modes 
in the reactive regime are those quasi-perpendicular 
to the beam. This is due to the large suppression 
caused by relativistic inertia along the longitudinal di- 
rection (|Fainberg et al.|[l~970h . However, as shown later, 
for quasi-perpendicular directions of the wave vector, the 
reactive regime is not applicable. 

When the approximation © is not valid, the growth 
rate is evaluated from a pole of the integrand in the dis- 
persion relation ([7]). namely 



Ik 



2ire 2 



df 

k ■ -± S(uj p - k ■ v) dp. (10) 



If, as is the case here, despite the energy spread 
the particles remain ultra-relativistic and, |v| = c, 
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Fig. 3. — Normalized growth rate as a function of kne/uip in- 
side the the resonant region, for values of k± where 7j. reaches its 
maximum values. 

ca n be assumed, the above i ntegral can be simplified 
to (jBreizman fc Mirnovlll970f ) 



n b (uj p \ 



dfx- 



(11) 



Here, the integration variable fj, is the angle between the 
particles and the beam direction, and 



p± = (u p /kc)m/k ± k ± /k^Jk 2 c 2 /oj 2 - 1), (12) 

(13) 



9(0) = ^- Pf(p,0)dp*= 

nb J 1 beam^-V 



e ao 2 



The second equality for g(9) in Eq. (fl~3|) is found to be 
a very good approximation based on results of the Monte 
Carlo model of the cascade. The integral for the growth 
rate in Eq. (fTTj) can be evaluated numerically. The qual- 
itative behavior of the growth ra te, 7fc, on the plane 
fc|l — k± is summarized in Fig. [2] ([Breizman fc Rvutov I 
|1971[ ) for a beam at a Gpc from the blazar. Outside the 
narrow resonant region of k-space, corresponding in the 
plot to the uniform color, the growth rate is effectively 
null. In the narrow resonant region around fcy = uj p /c, 
the growth can be positive (bright), negative (dark) and 
null, and for large enough values of k± , it carries the sign 
of, OJp/kn — c (see below). Within the resonant region, 
the growth rate as a function of fcii has typically two ex- 
trema, a maximum and a minimum. This is shown in 
Fig. [3] for values of k± of interest, i.e. where 7^. reaches 
its maximum values. The growth rate has its largest val- 
ues where k±c/uj p < 1, and decays rapidly in the opposite 
limit. This can be seen from Fig.|3]where 7^ is plotted for 
values of k± close to and much larger than oj p /c. Finally, 
we find that the growth rate, maximized with respect to 
fcii and as a function of k±, can be well approximated by 
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Fig. 4. — Normalized growth rate as a function of kuc/uip inside 
the the resonant region, for large values of k± where 74. starts to 
drop compared to its maximum value. 



the following expression by Breizma n fc Rvutov I (|1971[ ) 

lk -Up — -^79 . .9 , ,. 2 15 , (14) 



p Tn v A0 2 w2 + fc^c 2 ' 
which we will be using in the following. 

3.1. Fastest Growing Modes 

For particles of an ultra-relativistic beam with modest 
angular spread, A9 < 1, we can assume um ~ c and 
v± ~ cA8. If the energy spread, AE/E < 1, then the 
longitudinal velocity spread of the beam is 

AE 

Thus, for the angular spread and Lorentz factor char- 
acteristic of the blazars induced beam, the longitudinal 
velocity spread is negligible with respect to the perpen- 
dicular velocity spread. It turns out that, the first and 
second terms in Eq. (|15j) are comparable to within a fac- 
tor of a few, so for the sake of simplicity in the following 
we retain the second term only, neglecting any fudge fac- 
tor. We then combine Eq. flSJ, ([9]) and ([Poj) to obtain the 
condition for applicability of the reactive regime, i.e. 



2A9 < --± 

fen 



-i +4 

k 2 

h '\\ 



k 2 T 2 



kl 
k 2 



( 16 ) 

This threshold value for the angular spread of the beam 
is plotted in Fig. [5] as a function of the wave- vector ori- 
entation, for different values of the Lorentz factor of 
the beam, and the beam-to-IGM electron density ra- 
tio. For given beam characteristics, the nature of the 
instability depends on the orientation of the wave- vector 
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Fig. 5. — Threshold value of the angular spread for kinetic regime, 
as a function of the wave- vector orientation for the case of a beam- 
to-IGM electron density ratio of 10 -14 (red) and 10 -15 (cyan), 
respectively, and for a beam Lorentz factors T = 10 (solid) and 
5 X 10 5 (dash), respectively. 

with respect to the beam and modes of reactive and ki- 
netic character can appear simultaneously. In particular, 
modes close to perpendicular require a kinetic descrip- 
tion. This is because the beam velocity spread along the 
wave-vector cannot be neglected. 

Using again Eq. ([8|)- ([T"5"T) with the estimate of the beam 
angular spread and bulk Lorentz factor from the previ- 
ous Section, we find that the transition from reactive to 
kinetic regimes occurs for 



-i ~ xHT 5 
fc 



io- 15 
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(17) 



In the above estimate we have used a conservative value 
for the the ratio, nb/n v . However, this number may can 
be much smaller due to distance from the blazars and to 
redshift evolution of the EBL, so that virtually the reac- 
tive regime may never apply. In any case, for perturba- 
tion modes more aligned to the beam direction than the 
above value, the instability is reactive in nature, but the 
growth rate decreases due to the increasing relativistic 
inertia. In the opposite case, of larger k±/k, the insta- 
bility is kinetic and the growth rate also decreases be- 
cause of the increasing velocity spread along k. Thus, 
perturbations with fcj_/fc roughly as in Eq. (|17|) grow at 
the fastest rate, given by 
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(18) 



which can be actually derived either from the above re- 
active or kinetic expressions for the growth rate. 

4. BEAM STABILIZATION 

The generation of Langmuir waves induces particle- 
wave interactions which cause the particles to diffuse in 



momentu m space, eventu ally leading to the relaxation of 
the beam (|Melrose Ill989f ). In the very first phase of re- 
laxation, described by the quasilinear approximation, the 
generation of Langmuir waves is balanced by their drift at 
the group velocity, v g ~ Zvf /c, to regions of lower wave 
energy density. However, due to the smallness of the 
group velocity and to the large scale of the wave energy 
gradient, of order the pair beam distance to the blazar, 
the timescale associated to this process largely exceeds a 
Hubble time. Therefore in the quasilinear relaxation the 
growth of the Langmuir waves is uncompensated. 

4.1. Nonlinear Landau damping 

The main process that compensates the growth 

of Langmuir waves is induced scattering off plasma 

ions, also known as no n -linea r Landau damp- 
ing dTsvtovich fe Shapiro 



I197H iLesch fc Schlickeiser I 



119651 iBreizman et al. I 
19871) . In this process, a 
thermal ion, with characteristic velocity, v t i, interacts 
with the beat wave produced by two Langmuir oscilla- 
tions, w(fc), w(fc'), under the condition for Cherenkov 
interaction, i.e. 



ui(k) -w(jfc') = (fc - fc') 



(19) 



In a thermal plasma, the scattering process leads to 
damping o f the beat wave and a reduction of its 
frequency (| Kaplan fc Tsvtovich 1 119731 ). Due to their 
dispersion relation, Langmuir waves are consequently 
shifted towards regions of phase-space of lower wave- 
vectors. Eventually their phase speed exceeds the speed 
of light, cj/k > c, and they fall out of resonance with the 
beam particles. The change in w ave-vector suffered in 
one s cattering process is of order ()Kaplan fc Tsvtovich I 
[19731) 



Ak ~ uj„— = 50 x I 

Vie 



-Y 



V 10 3 K J 



(20) 



where v te is the electron thermal speed. Since for Lang- 
muir waves, fc ~ uip/c, the above equation implies 
Ak/k ^> 1, i.e. at the low temperatures characteristic 
of the IGM in voids, waves are kicked out of resonance 
in a single scattering e vent. The rate of induc ed scatter- 
ing off thermal ions is ([Breizman et al. Ifl97lh 



W 



(21) 



where W is the energy density of Langmuir waves. The 
latter evolves in time according to 



dW 
~dt 



= 2("frnax ~ Jnl)W. 



(22) 



Over time, a quasi steady-state condition will be estab- 
lished such that the growth of waves due to the beam 
instability is balanced on average by induced scattering. 
By equating the two competing rates we can obtain an 
estimate for the equilibrium value of W. Expressed in 
terms of the energy density of the beam, the latter is 



W, 



eq 



ni,Tm e c 2 
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3 x 10 5 
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10- 5 



The power associated with the emission of Langmuir 
waves by the beam is then dW/dt — 2W eq r y max - The 
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corresponding beam relaxation timescale is 
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(23) 



where redshift dependence is made explicit and 6 — 
—0.7 has been used. The above timescale is generally 
longer than the Hubble time, and even at close distances 
from the blazar, remains longer than the pairs cooling 
time on the Cosmic Microwave Background, tjc/c ~ 
xl0 6 (E ± /TeV)- 1 (l + z)- 4 yr. This suggests that the 
beam instability plays only a secondary role on the elec- 
tromagnetic shower and the beam dynamics and, as a 
consequence, the thermal history of the IGM and/or 
structure formation. 

4.2. Plasma Inhomogeneities 

In addition to induced scattering by thermal 
ions, plasma inhomogeneities can also induce 
loss o f resonance between particles and Langmuir 
wave s (jBreizman fc Rvutov 1119711 iNishikawa fc Rvutov I 
I1976T ) . The rate of change of the wave- vector due to 

inhomogeneities is given by dk/dt = -Vw. For Lang- 
muir waves the RHS of the above equation is mostly 
contributed to by density inhomogeneities. In the 
following we consider only the effects of inhomogeneities 
on the longitudinal component of the wave-vector. 
These are described by the equation 



dk\\ 
~~dt 



lUJp 

2 An 



(24) 



with Am = 7i v /(Vn v )|| , the length scale of the density 
gradient along the beam. 

Fig. |6] shows the range of scale lengths of IGM den- 
sity gradients as a function of IGM gas over-density, for 
three different cosmological redshifts, z = (top), 2=1 
(middle) and z — 3 (bottom). The data is extracted 
from a cosmological simulation of structure formation in- 
cluding hyd rodynamics, dark matter , and self-gravity as 
described in lMiniati fc Colellal (f2007h . A fiat ACDM uni- 
verse is adopted with the following cosmological param- 
eters: total mass density, normalized to the critical value 
for closure, f2 m = 0.2792, normalized baryonic mass den- 
sity, fit, = 0.0462, normalized vacuum energy density, 
Q A = 1 - TO = 0.7208, Hubble constant H = 70.1 km 
s _1 Mpc -1 , spectral index of primordial perturbation, 
n s = 0.96, and rms linear density fluctuation within a 
sphere with a comov ing radius of 8 fe" 1 M pc, erg = 0.817, 
where h = H /100 ([Komatsu et all 12009). The compu- 
tational box has a comoving size L = 50/i -1 Mpc, is 
discretized with 512 3 comoving cells, corresponding to a 
nominal spatial resolution of lOO/i -1 comoving kpc. The 
collisionless dark matter component is represented with 
512 3 particles with mass 6 x 10 5 /i _1 M Q . 

According to Fig. [51 the distance covered by the beam 
particles during the fastest growth time, ~ CHmax ~ 1 
kpc, is much shorter than the typical scale-length of den- 
sity gradients. This correspond to the case of regular, as 
opposed to random, inhomogeneities. In this case, the 
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Fig. 6. — Characteristic length scale of density gradient in the 
IGM as a function IGM gas over-density. The shaded area extends 
± one root-mean-squared value about the average. 

condition for wave excitation can be expressed as (Breiz- 
man & Ryutov 1971) 



Afcn 



In 



\dk»/dt\ 



> A c , 



(25) 



where A c is the Coulomb logarithm and Afeii the change 
in longitudinal component of the wave- vector allowed by 
the resonant condition (|6|). Eq. (|25[) states that, in order 
for the wave energy to grow, the resonance condition 
must be violated only after several growth times have 
lapsed. Using Eq. (|6]) (and neglecting the term AE/ET 2 ) 
to obtain, Afe|| < ^ A6 2 + k^A6, Eq. ([25]) can be solved 
to express the condition for wave excitation in terms of 
Ay, i.e. 



2uj p rib 



kuA6 



1.5xl0 7 kpcA c (i-') (^4? 

1 io 5 ) V io 15 



(26) 



where we have used, tu p = 16s _1 for the plasma frequency 
in voids. The above estimate gives an extremely stringent 
constraint on the minimum allowed scale-length of den- 
sity gradients required for the growth of Langmuir waves, 
which is orders of magnitude larger than the values re- 
ported in Fig. ([S]), even when using a very conservative 
value for the ratio, n v /nb- Therefore, such constraints 
cannot be possibly satisfied due to the presence of mat- 
ter fluctuations associated to cosmic structure, provid- 
ing another independent argument against the growth of 
Langmuir waves and the unstable behavior of the pair 
beam. 

5. CONCLUSION 
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We considered the stability properties of a low den- 
sity ultra relativistic pair beam produced in the inter- 
galactic medium by multi-TeV gamma-ray photons from 
blazars. The physical properties of the pair beam are 
determined through a Monte Carlo model of the electro- 
magnetic cascade. In summary we find that the com- 
bination of kinetic effects, non-linear Landau damping 
and density inhomogeneities appear to stabilize blazars 
induced ultra-relativistic beams over a Hubble timescale, 
so that the electromagnetic cascade remains unaffected 
by the beam instability. This implies that the lack of a 
bumpy feature at multi-GeV energies in the gamma-ray 
spectrum of distant blazars cannot be attributed such 



instabilities. Finally, heating of the IGM by pair beams 
appears completely negligible. 
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